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We study the structural transition from 'B' form of DNA to 'A' form of DNA using group theoretic 
methods. The transition is not of the order-disorder type and hence to construct a Landau kind 
of theory for the transition we define a higher symmetry and relevant order parameters. We also 
discuss the issue of all the conformations, observed experimentally during the course of transition, 
being fundamentally different or not. 



PACS numbers: 87.15.Zg 87.14.gk 

I. INTRODUCTION 

DNA can exist in many conformations, the most stable 
among them being 'A', 'B' and 'Z' forms. 'A', 'B' forms 
of DNA form right-handed double helix while 'Z' DNA 
forms left-handed double helix with varying number of 
turns per rotation. In this report we examine the transi- 
tion between 'B' and 'A' forms theoretically. 'B' DNA is 
the most commonly occurring DNA existing in highly 
humid conditions. It may be transformed to the 'A' 
form by methylation. The transition from 'B' form to 
'A' form was observed during the first x-ray studies by 
Franklin and Goslin[l] who characterized the change of 
the helix structure as ' a substantial re-arrangement of 
the molecule'. Experiments [2, 3] have been performed to 
observe the transition from 'B' form to 'A' form in a more 
detailed manner and this report is an attempt to explain 
the fundamental features of the experiments. Recent ex- 
periments^] to study the dynamics of the transition from 
'B' to 'A ' observe that the helical states are clearly sepa- 
rated from one another in terms of activation energy with 
the time scale of the dynamics of transition being of the 
order of ms. In the following paragraphs we will try to 
address the question of whether the transition from 'B' 
form to 'A' form is a phase transition as defined in statis- 
tical physics. The study of the 'B' to 'A' transition from 
the first principles by including the interaction between 
the various bases and their interaction with the environ- 
ment can be very challenging. The Landau analysis may 
provide a sophisticated method to study the phase tran- 
sition. 



II. SYMMETRY OF 'A' AND 'B' FORMS 

The Landau theory for the study of phase transition 
is considered as one of the important achievements of 



Landau. In Landau theory the free energy is written in 
terms of the order parameter which in turn depends on 
the symmetry of the phase having the higher symmetry 
through one of its irreducible representations. The free 
energy remains invariant under the symmetry operations 
of the higher symmetry. Phase transitions involving mul- 
tiferroics [5] complicated structures such as intercalation 
compounds [6] and antiferromagnets[7] have been studied 
using the Landau theory . In this report the Landau the- 
ory will be used to study the phase transition between 
'A' and 'B' forms of DNA. Experimental studies , among 
other techniques, use crystallography to study the struc- 
ture of DNA. The various conformations take on different 
crystallographic structures[2,3] corresponding to different 
space groups. The space groups depend on the detailed 
symmetry of the molecule, length of the sequence etc. 
The use of space groups may not be the most reliable 
method to study the transition theoretically After all, 
the DNA exists as a single macro molecule and we would 
prefer to treat it like that. 

In the present approach the DNA is considered to be 
double helix with Watson-Crick bonding between the two 
strands. We would not go into the detailed coordinates 
of the atoms in the two forms but would like to explain 
the transition from the point of view of change in symme- 
try resulting in the difference between 'A' and 'B' forms. 
With this objective in mind, we consider the DNA macro 
molecule to be infinitely long along the axis. Hence , it 
can be looked upon as a transition where the line group 
changes [8]. The view along the axis of 'A' and 'B' DNA is 
shown in Figl[9]. Apart from a pure translational sym- 
metry 'a' along the axis, 'B' form has the screw symme- 
try 10i. The rotation part of the group is cyclic with 
order 10(Cio)[9]. 'A' form of DNA has screw symmetry 
Hi. The rotation part of the group is cyclic with order 
ll(Cn). In the subsequent analysis we employ group 
theoretic arguments to construct the order parameter in 
order to write down the free energy. The irreducible rep- 
resentation for the rotation part, say of order a, may 
be deduced from the relation z a — 1 as the group is 
cyclic. So z — e 2- ^ where n = 0, l...,a — 1. Since the 
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remains invariant under a subgroup of transformation in 
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the ordered phase. Hence ^^f + hz = 2-kA , A being an 



integer. It leads to 



k 



2tt(10A - n) 



(2) 



One may repeat the same arguments for 'A' DNA to 
obtain, 



A; = 



2tt(11A - n) 



(3) 



FIG. 1: Axial view of 'A' and 'B' forms 

groups(Cio,Cn) are abelian, all the representations are 
one dimensional . D n represents the n th irreducible rep- 
resentation. 

The basis for the representations is given by e~ m< ^ 
where n = 0,1, ..9 for Cio and n — 0,1, 2. .10 for C\\. 
R(Acp) be an operator that rotates the coordinate sys- 
tem by Ac/). U = % (say), then i?(A0) e - m * = 
eTe - ™^, n being the index of representation. Thus 
e^S 2 " is a representation ofrotation by 27r/10 with index 
n. We can thus calculate the representations of C\q and 
C\\ corresponding to the various irreducible representa- 
tions. The irreducible representation of the line group is 



given by e 



e being the basis for the represen- 



tation of translation along the DNA axis. The represen- 
tation is labeled by the pair (n, k). 



III. CONSTRUCTION OF THE ORDER 
PARAMETER 

The order parameter, which is used to describe order- 
disorder transition, is the expectation value of an op- 
erator which does not remain invariant under the full 
symmetry group, taking the value in the disordered 
phase and a finite value in the ordered phase. It, though, 
remains invariant under a subgroup of transformations. 
In order to introduce the order parameters in the prob- 
lem we need to define disordered phase ( high symme- 
try phase). We consider 10i('B' DNA) and lli('A' 
DNA) as the ordered phases in comparison to the dis- 
ordered phase(high symmetry phase) given by the sym- 
metry (10 x 11 x m)i, m being some integer. This is so 
as they are subgroups of the high symmetry phase. We 
take m = 1 in this case for the sake of convenience. Let 
) T ( i§) = °- 9 ) be a symmetry operation of 'B' 
DNA with R standing for the rotation and T standing 
for the translation. 
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We want e i° e u to be unity since the order param- 

ptpr transforms as sranp irrprlnri hip rpnrpspntation and 
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FIG. 2: Schematic diagram representing the Landau 
theory for transition between 'A' and 'B' forms of 
DNA 

We note that n = 0, ..109 since the higher symmetry 
is denoted by (10 x We take n = 1, A = 1 for 

convenience. For 'B' DNA k = — and for 'A' form 

a 

k = though there could be different values of k cor- 
responding to different values of n and A. We would 
work in the space of basis (fci,ni) and (£2,^2) where 
fei = ^ ,m = 1 and k 2 = ^,n 2 = 1. There is no 
mixing between (Jfci,7ij.) and (fc2,?i2) > so the system is 
defined by two order parameters a a and otb transform- 
ing respectively as the representations determined by the 
above basis. (a Q ,0) is invariant under the symmetry op- 
erations of Hi — <^a remains invariant due to (3) and 
remains so under the operation — thereby representing 
the 'A' DNA phase, phase. For similar reasons (0, a ) 
being invariant under the operations of 10i represents the 
'B' DNA phase . (0,0) being invariant under the larger 
symmetry group represents the melt or the higher sym- 
metry The picture of the phase transition is presented 
in Fig2. We try to establish the transition between B 
DNA and A DNA by looking at the transitions, possibly 
fictitious, between the disordered phase and the ordered 
phases. We consider a a and at to be complex quanti- 
ties. Having constructed the two order parameters we 
may now attempt to write down the free energy of the 
system. It may be seen that under a general operation of 
(10 x ll)i a a and a& picks up a phase. Hence a very gen- 
eral form of the free energy in the form of a polynomial 

in rv axiH rv* until the miartic orHpr. invariant under the 
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operation of the higher symmetry group, may be written 
down. 
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a b 



- -<?(! «« I 2 



«6 I 2 ) 



+Ui | a a | 4 +u 2 | a b 4 +2w 22 a a | 2 | a b 



(4) 
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. It may be seen here that the coeffi- 
a b | 2 have been split among r and g as 



where Wi,it2 
cients of | a a 

they separately carry different information. The transi- 
tion from 'B' DNA to 'A' DNA is carried out through the 
addition of a chemical agent (c). In the above equation 
we put r oc (T — T*) and g oc (c— c*), T* and c* (chemical 
cone) are some constants and T, not necessarily the tem- 
perature, may not be a realizable parameter. M22 is the 
term that is responsible for the coexistence of phases. 
The reason for this would be clear when we discuss the 
solutions to (4). The minimization of the free energy 
gives the following combinations for the possible values 
of a a and a b . 



1=0,1 oi, |=0 



a a = I, I a b \yt 



a a |^ 0, 1 at |= 



a a |^0,| a b |^0 



(5-a) 



(5-b) 



(5-c) 
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We first consider the case where U22 > 0. A complete 
analysis involving the stability of the solutions in the 
r — g plane is shown below[10]. Fig3 corresponds to the 
case uiu 2 — u 2 2 2 < and Fig4 corresponds to the case 
U1U2 — U22 2 > The dotted lines in the figures denote 
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FIG. 3: Phase diagram corresponding to the free 

energy given by eqn(4) in the r — g plane for 
U1U2 — U22 2 < 0. This diagram corresponds to the 
transition between 'A' and 'B' forms. 

second order transition while the wia'P'lv line denotes first 
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FIG. 



Phase diagram corresponding to the free 
— g plane for 



energy given by eqn(4) in the r 

U1U2 — U22 2 > 0. 



order transition line. Being on the wiggly line signifies 
the state of coexistence as shown in Fig3. In Fig3 we 
observe that the disordered phase (| a a = |, | a b = |) 
is separated from the ordered phase('B' form and 'A' 
form) by the second order lines while the ordered phases 
themselves are separated by first order line. The sec- 
ond order and the first order lines meet at the bi critical 
point. In Fig4 we observe the emergence of another or- 
dered phase(| a a 1 7^ 0, 1 a b \fi 0) separated from the 
other two ordered phases by second order lines. If m 2 2 
is sufficiently positive so that the condition for Fig3 is 
satisfied, we have coexistence, otherwise not. Consider- 
ing U22 —> 0, the condition for Fig4 is easily satisfied and 
hence no phase coexistence as should be the case. For 
Mi > 0,1*2 > and 1*22 > we get consistent solutions 
for (4). But if U22 < 0, it can be shown that the re- 
quirement of the stability of the solutions gives rise to 
inconsistent picture. This result is entirely an analyti- 
cal outcome of (4) and is independent of DNA. We have 
used group theoretic arguments to construct the order 
parameters and write the free energy expression. The 
free energy gives certain phase diagrams which can be 
explained in the light of experimental results. It may be 
noted here that any transition taking us from the group 
Pi to gi may be represented by a similar theory. 



IV. INTERPRETATION OF THE ORDER 
PARAMETER 



Having constructed the order parameters we would 
now like to identify them. In condensed matter sys- 
tems, we may have multicomponent order parameter de- 
termined by the possible irreducible representation of the 
point group of the relevant vector k in the reciprocal 
lattice a.s well as the star of k. . Tn the case of antifer- 
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romagnets[6] we may , for example, have the magnetiza- 
tion perpendicular to the antiferromagnetically coupled 
ferromagnetic planes or situated in the plane. For the in- 
tercalation compounds, the components of the order pa- 
rameter are in a way related to the Fourier components 
of the density of the alkali atoms corresponding to the 
stars of k . In the present case, however, we are dealing 
with two single component order parameters for a sys- 
tem whose structure does not fall under the category of 
Bravais lattice. Hence the above prescription would not 
be followed while identifying the order parameter. Being 
a structural phase transition, we would relate the order 
parameters to densities. For the double helix DNA, Hi, 
10i or (lOx ll)i symmetry uniquely defines the sites on 
the strands. p a denotes the densities of sites satisfying 
Hi symmetry and pb the density of sites satisfying 10 1. 
We identify the modulus of the order parameters as re- 
lated to p a — p and pb — pb = if p a > pb and vice versa. 
In the 'A' phase pb = and in the 'B' phase p a = 0. 
In the case of p a = Pb (disordered phase), we have for 
both the order parameters p a — Pb = 0. With the above 
definition of the order parameters, we are able to identify 
the two ordered phases and the the disordered phase, as 
depicted in Fig3. 

Experiments that have been carried out to study the 
phase transition [2,3] report conformations in which both 
'A' type and 'B' type features are there apart from only 
'A' type or 'B' type features. We interpret this as coex- 
istence of phases. This may suggest that the transforma- 
tion is hrst-order in nature as depicted in Fig3, coexis- 
tence being a characteristic of first order transition. The 
two phases coexist for the order parameters jump dis- 



continuously at the first-order line. In first order phase 
transition when one goes from one phase to another there 
is jump in the order parameter at the transition point . 
Similarly when one goes back to the same phase there is 
a jump. This only implies coexistence of phases at the 
transition point. Further, it may be noted that the DNA 
molecule, under experimentation in [2,3], being not a 
long molecule, has a finite region of coexistence in the r-g 
plane whereas Fig3 shows a line(wiggle) representing the 
coexistence valid for infinitely long molecule. | a a \ j£ , 
I ctb \fi is the most ordered phase and in our problem 
should reflect a phase that has the lowest symmetry.But 
such a phase which will have a lower symmetry than both 
10i and Hi is not possible. This perhaps tells us that 
the Fig3 represents the likely phase picture and Fig4 may 
be ruled out. 

If we refer to the experimental results of [2,3], it is men- 
tioned that the transition from 'B' DNA to 'A' DNA 
takes place through 14 conformations. Further, the ex- 
perimental data in [2,3], where they have measured cer- 
tain quantities like tilt angle, axial distance etc do not 
actually tell us whether there has been a phase transi- 
tion or not. But in the light of the above analysis, it 
may be concluded that the transition from 'B' to 'A' is 
a first order phase transition with all the conformations 
being not fundamentally different but belong to different 
phases('B' and 'A') and the region of coexistence. 
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